Short-range magnetic correlations in Tb5Ge4 
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We present a single crystal neutron diffraction study of the magnetic short-range correlations 
in Tb5Ge4 which orders antiferromagnetically below the Neel temperature Tat ~ 92 K. Strong 
diffuse scattering arising from magnetic short-range correlations was observed in wide temperature 
ranges both below and above Tm- The antiferromagnetic ordering in Tb5Ge4 can be described 
as strongly coupled ferromagnetic block layers in the ac-plane that stack along the 6-axis with 
weak antiferromagnetic inter-layer coupling. Diffuse scattering was observed along both a* and 6* 
directions indicating three-dimensional short-range correlations. Moreover, the g-dependence of the 
diffuse scattering is Squared-Lorentzian in form suggesting a strongly clustered magnetic state that 
may be related to the proposed Griffiths- like phase in Gd5Ge4. 

PACS numbers: valid numbers to be inserted here 



I. INTRODUCTION 

Tb5Ge4 and Gd5Ge4 belong to the rare earth 
R5(Sia;Gei_2;)4 series compounds. These materials 
exhibit large magnetocaloric (MC) efFects^i^ii^i^ and 
are currently attracting much attention for their 
potential appUcation in magnetic refrigeration^i^i^ii. 
Both Tb5Ge4 and Gd5Ge4 are rich in magnetic 
propertie a^'^'^^i^^'^^i-'^^'^'^i^^'"'^^ and are believed to play a 
key role in understanding the underlying physics of the 
R5(Sia;Gei_2;)4 systems. They both crystallize in the 
Sm5Ge4-type crystallographic structure and adopt the 
same magnetic space group, Pnm' a^^^^'^- ' ^^ . Tb5Ge4 and 
Gd5Ge4 undergo long range antiferromagnetic (AFM) 
transitions at ~ 92 K and ~ 127 K, respectively. The 
magnetic structure of Tb5Ge4 and Gd5Ge4 consist of 
Tb/Gd-rich block layers in the ac-plane that stack along 
the 5-axis with strong ferromagnetic intralayer and weak 
AFM interlayer interactions. Gd5Ge4 has a coUincar 
AFM structure with the magnetic moments lying within 
the block layers along the c-axis. Tb5Ge4 orders in a sim- 
ilar fashion, although the single-ion anisotropy results in 
significant canting of the moments at low temperature. 

Figure [1] shows the magnetic susceptibility of Tb5Ge4 
along all three crystallographic axes measured with a 
Quantum Design SQUID Magnetic Properties Measure- 
ment System. Two phase transitions were observed: the 
AFM transition at T/v 92 K and a second phase tran- 
sition at ~ 55 K. The ~ 55 K transition was attributed 
to a spin reorientation transition similar to that pro- 
posed for Gd5Ge4 at ~ 75 K— . It has been suggested 
that the spin reorientation transition in both Gd5Ge4 
and Tb5Ge4 may arise from the delicate competition be- 
tween the magnetic anisotropy from the spin-orbit cou- 
pling of the conduction electrons and the dipolar inter- 
action anisotropy In addition to the AFM and spin 
reorientation transitions, it is of particular interest that 
significant magnetic short-range correlations (SRC) was 
suggested in Gd5Ge4 at temperatures both below and 
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FIG. 1: (Color online) Temperature dependencies of the mag- 
netic susceptibility measured with applied magnetic field (H 
= 1000 Oe, zero field cooled) parallel to all three crystallo- 
graphic axes. 



above T/v based upon the low magnetic field dc magneti- 
zation and ac magnetic susceptibility measurements^^. It 
has been interpreted as evidence of a Griffiths-like phase 
similar to the one observed in Tb5Si2Ge2, investigated 
by small-angle neutron scatteringi^. A Griffiths phase 
(GP)S2. is a nanoscale magnetic clustering phenomenon 
that is usually associated with competing magnetic inter- 
actions in the system21i22iiS^. It is possible for a Griffiths- 
like phase to exist in Gd5Ge4 due to the competition 
between FM and AFM interactions in Gd5Ge4, FM in- 
teractions within the layers, and either AFM or FM in- 
teractions between the layers (small structural changes 
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FIG. 2: (Color online) Longitudinal scans (in LoglO scale) along the (OfcO) direction measured at several temperatures (a) 
below Tjv, (b) above Tjv, indicate strong AFM SRC with broad magnetic scattering observed around k — odd. The remaining 
scattering intensity of (010) and (030) at 200 K is likely due to multiple scattering. 



or applying relatively low magnetic field comparing to 
T/v can switch the interlayer interactions between AFM 
& FM interactions therefore switch the low temperature 
phase between AFM order or FM order) . However, stud- 
ies of Gd5Ge4 have been hampered by the large neutron 
absorption cross-section of gadolinium. Tb5Ge4 exhibits 
similar magnetic properties as to Gd5Ge4 hence it is an 
ideal candidate for neutron scattering studies of magnetic 
SRC displayed in these compounds. We report here neu- 
tron diffraction studies of the AFM phase transition and 
the magnetic SRC in Tb5Ge4. 



II. EXPERIMENTAL DETAILS 

A large Tb5Ge4 single crystal (~ 5 grams) was used 
for the neutron diffraction experiment. The single crys- 
tal was grown at the Materials Preparation Center— us- 
ing the Bridgman technique as described in Ref. [2^. 
The mosaic of the crystal is 0.41(3)° along the a-axis 
and 0.63(2)° along the 6-axis as determined by the full 
width at half maximum of the (400) and (060) Bragg 
peak rocking curves. The crystal was mounted on a thin 
aluminum post, oriented in the (hkO) scattering plane, 
and sealed in a helium filled aluminum sample can. A 
closed-cycle Helium refrigerator (Displex) was used which 
allows accurate temperature control between 10 K and 
300 K. The experiments were performed using the HBIA 
triple-axis spectrometer located at the High Flux Iso- 
tope Reactor (HFIR) at the Oak Ridge National Labo- 
ratory (ORNL). The HBIA spectrometer operates with 
a fixed incident energy, Ei = 14.6 meV using a double 
pyrolitic graphite (PC) monochromator system. Two 
highly oriented PC filters were mounted before and after 



the second monochromator to significantly reduce higher 
order contaminations of the incident beam (i.e., 1^/2 = 
IO^^'/a). a coUiniation of open-40'-sample-'i0'-68' was 
used throughout the experiment. All data have been nor- 
malized to the beam monitor count. 



III. RESULTS AND DISCUSSIONS 

Figure 2 compares longitudinal scans measured along 
the (OfcO) direction at selected temperatures both be- 
low (Fig. m (a)) and above (Fig. [2] (b)) Tn- Strong 
AFM magnetic reflections with k — odd integer were ob- 
served below T^v consistent with the magnetic structure 
of Tb5Ge4. At low temperatures, the (010) and (030) 
magnetic reflections are superimposed on weak and very 
broad Lorentzian-shaped diffuse scattering peaks aris- 
ing from magnetic fluctuations. Below T^, the diffuse 
scattering increases with increasing temperatures as il- 
lustrated in Fig. [5] (a). The diffuse scattering is strongest 
at T^r and then weakens as the temperature increases 
above Tn as shown in Fig. [2] (b) as expected for typical 
critical behavior. The data indicate strong diffuse scat- 
tering around the strong magnetic reffections over a wide 
temperature range. 

We characterized the Tjv 92 K AFM transition with 
an order parameter measurement. Fig. [3] depicts the or- 
der parameter of Tb5Ge4 as measured by monitoring the 
strong magnetic reffection (030) as a function of temper- 
ature. The integrated intensity was obtained by fitting 
the (030) rocking curve measured at each temperature 
to a Lorentzian function with a constant background. A 
fit of the order parameter to a power-law I(T)=Iq[(Tn- 
V/T^Jf^ yields Tn « 91.38 ± 0.05 K and /3 « 0.20 
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FIG. 3: (Color online) Tb5Ge4 order parameter, integrated 
intensity of the (030) magnetic peak as a function of temper- 
ature. The red line is the fitting of the order parameter data 
to the power law as described in the text which yields Tjv 
^ 91.38 ± 0.05 K and /3 0.20 ± 0.01, and the dashed and 
dotted lines are calculations using the same fitting parame- 
ters but simply replacing the 13 value to be the values of a 2D 
Ising (/3=0.125) and a 3D Ising (/3=0.326) system. 
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FIG. 4: (Color online) Diffuse scattering intensity (88 K - 8 
K) vs. k and h around (030) constructed from a series of grid 
scans measured at 8 K and 88 K. The 8 K data is subtracted 
from the 88 K data to eliminate contributions from (030) 
magnetic Bragg reflection. 



± 0.01, where (5 is the critical exponent. The fitting re- 
sult is plotted in Fig. [31 solid red line, in comparison 
to calculations using the same fitting parameters but re- 
placing 13 with the values of a two-dimensional (2D) Ising 
(/3=0.125) and a three-dimensional (3D) Ising (/3=0.326) 
system2&. The obtained critical temperature Tjv is in 
good agreement with the magnetic susceptibility result. 
The yielded (3 value ^ 0.20 is between the theoretical val- 
ues of a 2D and a 3D Ising system, which suggests that 
the dimensionality of Tb5Ge4 is intermediate between 2D 
and 3D consistent with its layered structure. 

The longitudinal scans shown in Fig. [2] indicate strong 
diffuse scattering along the 6-axis. In order to see how 
the diffuse scattering is distributed in the (hkO) plane, a 
series of grid scans around (030) were performed at two 
temperatures, 8 K and 88 K. The 8 K data were sub- 
tracted from the 88 K data to eliminate contributions 
from the (030) magnetic Bragg reflection. Fig. |4] is the 
contour plot of the subtracted diffuse scattering inten- 
sity vs. h and k. It shows that the diffuse scattering 
extends along both h and k directions. The diffuse scat- 
tering intensity is strong around (~0.1 ^^3.1 0) and its 
equivalent positions. No strong anisotropy is observed 
suggesting the magnetic correlations associated to the 
diffuse scattering is not restricted to the FM block layer 
(ac-plane) but is rather more three-dimensional exhibit- 
ing AFM critical scattering behavior. 

The diffuse scattering was studied in detail as a func- 
tion of temperature. Wide transverse scans along the 
h direction were performed at (0 1.12 0) to reduce con- 
tributions from the nearby (010) strong magnetic reflec- 
tion. At 10 K, the FWHM (full width at half maximum) 
peak width of the (010) is about 0.09 [r.l.u] (reciprocal 



lattice unit). Therefore we may attribute the scatter- 
ing intensity measured at (0 1.12 0) to diffuse scatter- 
ing from magnetic correlations. Fig. [5] (a) shows typical 
scans at different temperatures. The scattering intensity 
of the strong diffuse scattering peak first increases with 
increasing temperature up to T/v and then decreases with 
further increasing temperature above T/v. In general, the 
g-dependence (here q — h or k) of diffuse scattering can 
be well described by the following equatio n^-i -- - , a sum of 
a Lorentzian function and a Squared-Lorentzian function 
plus constant background (BG) , 
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The first Lorentzian term is the conventional critical 
scattering component representing an Ornstein-Zernike 
form, i.e. exp~'""/r, for the magnetic correlations^^. The 
second Squared-Lorentzian term^S is generally consid- 
ered to arise from static or frozen spin clusters within 
which the spin correlations decrease more gradually as 
exp"'"'. Surprisingly, by fitting the data to a Lorentzian 
-I- BG function only, a Squared-Lorentzian -|- BG function 
only, and the sum of Lorentzian plus Squared-Lorentzian 
function as described in Eq. [1] respectively, we found 
that the diffuse scattering data can be best described 
by the Squared-Lorentzian function only plus constant 
background. Adding an extra Lorentzian term docs not 
improve the quality of the least-squares fit to data. The 
comparisons between the fits to a Lorentzian only and a 
Squared-Lorentzian only are shown in Fig. [51 Note that 
the diffuse scattering peak is very broad, thus the resolu- 
tion effect can be neglected and no resolution corrections 
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FIG. 5: (Color online) Diffuse scattering along the /i-direction measured at (0 f .12 0). (a) Representative scans measured at 
different temperatures. The solid curves are least-squares fits to a Squared-Lorentzian function of h as described in the text; (b) 
Comparison of the fits of the 95 K data to a Lorentzian + BG (blue dashed line) function only and a Squared-Lorentzian -f- BG 
(red solid line) function only; (c) vs. T, where the normalized is obtained from least-squares fits to a Lorentzian + BG 
only and a Squared-Lorentzian -I- BG function only; (d) Integrated intensity vs. temperature; (e) Peak width vs. temperature; 
The dashed black lines in (d) and (e) are fits of the T > Tjv data to a power law as described in the text. 



are applied in the data analysis. Fig. [5] (b) compares the 
fits of the 95 K data to Lorentzian -I- BG (blue dashed 
fine) only and Squared-Lorentzian -I- BG (red solid line) 
only. This indicates that the 95 K data can not be de- 
scribed by a Lorentzian -I- BG function. On the other 
hand, as shown in Fig. [5] (a), the solid curves are fits to a 
Squared-Lorentzian + BG only that adequately describe 
the data for the measured g-range at all temperatures. 
The comparison of the obtained normalized from these 
two fittings is shown in Fig. [5] (c). It clearly shows that 
the data are better captured by a Squared-Lorentzian 
form, the Lorentzian line shape does not give a good de- 
scription of the data, particularly at temperatures near 
Tjv as illustrated in Fig. [H] (b). 

The integrated intensity and the FWHM obtained 
from least-squares fits to the data with a Squared- 
Lorentzian function are plotted in Fig. [S] (d) and (e). 
Two features are observed in the integrated intensity data 
(Fig. [5] (d)). The small kink at ~ 55 K is associated with 
the spin reorientation transition, and the peak at ~ 92 
K is associated with the AFM transition. Both temper- 
atures agree well with the magnetic susceptibility data. 



Despite the small kink at ~ 55 K, the ~ 92 K peak is 
nearly symmetric indicating strong critical fluctuations 
at T/v that die off as one moves away from Tn in cither 
direction. Above T^r, the integrated intensity data can 
be fit to a power law I(T)=Iq[(T-Tn )/Tn yielding 
Tn « 89.94 ± 1 K and /?' « 0.22 ± 0.02 (dashed line in 
Fig. El (d)). The obtained /3' value agrees to the (i value 
obtained from the fit to the order parameter data. As il- 
lustrated in Fig. [5] (e), the correlation length of the SRC 
remains relatively constant below T/v as indicated by a 
nearly constant peak width. Above T/v, the peak width 
increases as expected as the correlation length decreases 
with increasing temperature. A fit to the T > T/v peak 
width data to a power law S,(T)=io[(T-TN )/Tn )J" gives 
T/v ~ 91.7 ± 1 K and i; 0.77 ± 0.06 (dashed line in in 
Fig. [5] (e)), the v value is between the theoretical values 
of a 3D Ising [v = 0.6312) and a 2D Ising {v = 1)^ con- 
sistent with the order parameter measurement results. 
The Squared-Lorentzian peak widths obtained from fits 
to (010) longitudinal and transverse scans are also shown 
in Fig.[5](e). The Lorentzian-squared lineshape provides 
the best fit along both the h and k directions, indicating 
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that the correlations in the spin clusters extend both in 
the block layers and between the blocks. 

Our neutron diffraction study reveals strong diffuse 
scattering in Tb5Ge4 that persists to temperatures well 
above Tm- A detailed study of the peak shape indicates 
it is not conventional critical scattering with a Lorcntzian 
shape but shows a Squared-Lorentzian peak shape. As 
described in Ref. 29 and Ref. 30, the Squared-Lorentzian 
term arises if the pair correlation function falls off as 
exp"'''", which is characteristic of a spin-cluster state. 
Although the diffuse scattering is Squared-Lorentzian in 
form providing evidence of a clustered magnetic state in 
Tb5Ge4, we believe that the diffuse scattering observed in 
Tb5Ge4 is quite different from the proposed FM GrifRths- 
like phase in Gd5Ge4 (inferred from dc/ac magnetization 
and magnetic suscpetibility studies) for the following 
reasons. (1) As depicted in Fig.[2l at temperatures both 
below and above Tjv, the diffuse scattering is peaked at 
odd values of k (AFM wavevector) only, indicating that 
it is associated with AFM fluctuations. (2) The inte- 
grated intensity of the diffuse scattering also behaves like 
that typical of magnetic critical fluctuations, with a di- 
vergence of the correlation length at Tm- (3) The critical 
exponents obtained by fits of the diffuse scattering in- 
tegrated intensity and peak width to a power law are 
consistent with the values obtained from AFM order pa- 
rameter measurements. (4) Quasi-elastic measurements 
indicate the diffuse scattering is static in origin. Our neu- 
tron diffraction data indicate that the diffuse scattering 
observed in Tb5Ge4 exhibits behaviors of AFM critical 
fluctuations despite the Squared-Lorentzian peak shape. 



The fact that the peak shape of the diffuse scattering 
is not a Lorcntzian, as expected for normal critical 
scattering, is interesting and should not be left without 
a discussion. Here we consider two possibilities that 
may affect the diffuse scattering peak shape. (1) The 
Squared-Lorentzian peak shape may be intrinsic, i.e. 
related to the Griffiths-like phase, formation of which 
has been discussed in Refs. [l^ andfl^ (2) The unusual 
peak shape may also arise from some extrinsic effects, 
for example impurities in Tb5Ge4. It has been reported 
that R5(Sia;Gei_j,)3-type impurity phases, seen as 
very thin plates that are scattered through the bulk 
of R5(Sia;Gei_2,)4 samples, are present in all studied 
compounds of this series regardless of R^. Our data 
show that TbsGcs impurity phase is also present in the 
studied Tb5Ge4 crystal. It is possible that the Tb mag- 
netic sublattices of Tb5Ge4 are disrupted by the TbsGea 
impurities resulting in spin-clusters in Tb5Ge4 which 
give rise to the Squared-Lorentzian diffuse scattering 
peak shape. 
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